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Abstract 

We study weakly coupled SU{N) A/" = 4 super Yang-Mills theory on M x S"^ at infi- 
nite A^, which has interesting thermodynamics, including a Hagedorn transition, even at 
zero Yang-Mills coupling. We calculate the exact one-loop partition function below the 
Hagedorn temperature. Our calculation employs the representation of the one-loop dilata- 
tion operator as a spin chain Hamiltonian acting on neighboring sites and a generalization 
of Polya's counting of 'necklaces' (gauge-invariant operators) to include necklaces with 
a 'pendant' (an operator which acts on neighboring beads). We find that the one-loop 
correction to the Hagedorn temperature is SlnTu = +X/8n'^. 
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1. Introduction 

The past several years have witnessed a tremendous amount of effort invested into 
the careful study of A/" = 4 super Yang-Mills (SYM) theory at large from a number of 
complementary approaches. One motivation for much of this work is the fact that this 
theory is believed to provide, via the AdS/CFT correspondence, the simplest context in 
which we might hope to understand how to solve large A^ gauge theories in four dimensions. 
Optimistically, the apparent integrability of the string theory dual to A/" = 4 SYM theory 
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suggests that it might be possible to calculate all physical quantities (at least at infinite 
A^) as exact functions of the 't Hooft parameter A = Qym-^- The successful calculation of 
the circular BPS Wilson loop for all A provided an early realization of this hope. 

More recently, significant progress has been made on the problem of calculating anoma- 
lous dimensions of gauge theory operators. Two related approaches relevant to this problem 
include the study of semiclassical string solutions following [Q, and the study of the di- 
latation operator directly as the Hamiltonian of an integrable spin chain following H,^. 
Comprehensive reviews of the most recent progress on these approaches can be found in 
Ij^,^ and 0. BMN operators provide a prime example of a class of operators whose anoma- 
lous dimensions can be calculated for all A [|10|,|Tl[]. Very recent work with partial results 
on the problem of summing up all orders in A includes [p!2| , p!3[| . 

In this paper we continue in the fine tradition of chipping away at A/" = 4 SYM theory 
from a variety of angles. Our present interest lies in the partition function of the theory 



on M X 5" at infinite A^, which displays interesting thermodynamic behavior |14], including 



a Hagedorn transition [|T5[, even at zero 't Hooft parameter. The partition function can 
be calculated exactly at A = by simply counting gauge-invariant operators using Polya 
theory (see for example ||T5| , |T6| , [T7| , |T8| ) . Here we calculate the one-loop correction to the 



partition function below the Hagedorn temperature and find the result 



, oo 



loop ^ ^y]-2 



1 - z{Xk) 



(1.1) 



where Xn = c<;"'"'"^e~^"', uj = e^'^* (so that cu™"/^ = ±1 depending on whether m is even or 
odd), and the function z{x) is given below in (2.4). The remaining quantities {D2{x)) and 
{PD2{w,y)) are traces of the one-loop dilatation operator D2 acting on two neighboring 
fields inside an operator. We evaluate these traces explicitly in (5.14) and (5.17) below. We 
also obtain from ( |1 . Ij ) the one-loop correction to the value of the Hagedorn temperature, 

Th(A) =Th(0) (^l + ^ + ---^ , Th(0)^ 0.380 (1.2) 

(measured in units where the radius of the is one), which is presumably another example 
of an interesting physical quantity which we might hope to one day calculate as an exact 
function of A. 

We begin in section 2 by reviewing how the free partition function may be calculated 
by first assembling the partition function for the elementary fields into a partition function 
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for single-trace operators, and then into the fuU partition function for multi-trace operators. 
We discuss the general structure of the one-loop correction to the partition function and 
comment on the role of operator mixing and effects. In section 3 we rephrase the 
problem of counting gauge theory operators into the language of spin chains and introduce 
a generalization of Polya's necklace problem to what we call necklaces with a 'pendant'. 
These are necklaces with some local operator O, such as a spin chain Hamiltonian, inserted 
at neighboring beads on the necklace. We derive a general formula expressing a partition 
function for such necklaces in terms of some basic quantities {0{x)) and {PO{w, y)) which 
are easily obtained from O. In section 4 we apply this machinery to various familiar 
subsectors iSU{2), SO{6), and SL{2)) of the TV = 4 SYM theory. This section is mostly 
a warm-up for section 5, where we calculate {D2{x)) and {PD2{w,y)) for the one-loop 
dilatation operator D2 of the full A/" = 4 SYM theory. Finally in section 6 we present our 
final results for the one-loop corrections to the single- and multi-trace partition functions 
of AT = 4 SYM. 



2. The Af = 4 SYM Partition Function 

We begin in this section with a discussion of weakly-coupled AT = 4 SYM theory on 
M X S""^ following [1T5|JT^] . Since there is more than one way to calculate the free partition 
function, we review here that derivation which best sets the stage for our calculation of 
the one-loop correction in the following sections. 

2.1. Initial considerations 

In general, the partition function is given by 

Z{P)^Tr[e-f'% (2.1) 

where P is the inverse temperature and H is the Hamiltonian of the theory on M x S"^. It 
is convenient to introduce the bookkeeping parameter 

x = e-^ = e-i/T, (2.2) 

which ranges from x = (zero temperature) to x = 1 (infinite temperature). According to 
the state-operator map there is a one-to-one correspondence between states of the theory 
on M X iS"^ and gauge-invariant operators on R^. The Hamiltonian H onM-X is identified 
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with the dilatation operator D on the plane, so counting states weighted by x to the power 
of their energy is equivalent to counting gauge-invariant operators weighted by x to the 
power of their dimension on the plane. Therefore the partition function can be written as 

Z{x) = Tr[x% (2.3) 

where we implicitly set the radius R of the to one. The dimensions of quantities such 
as energy and temperature can be restored by affixing the appropriate power of R. 

We start by calculating Z{x) at tree level, so the dilatation operator reduces to D = 
-Do 5 which just counts the engineering dimension of an operator. To calculate ( |2.3|) all we 
need to do is write down a complete basis of operators and count them. The most general 
gauge-invariant operator can be written as a linear combination of operators with a definite 
number of traces, and the most general /c-trace operator can be expressed as a product of 
k single-trace operators (keeping in mind that separate traces behave as identical particles 
and are subject to the appropriate Bose or Fermi statistics). Therefore, a complete basis 
for arbitrary gauge-invariant operators follows naturally after we specify a complete basis 
of single-trace operators. 

2.2. TheU = A alphabet 

A single-trace operator is a product of the fields {(j)^ , A", F^i,) of A/" = 4 SYM theory 
and their covariant derivatives. Covariant derivatives must always be symmetrized with 
traces removed, since antisymmetric derivatives can be replaced by field strengths and 
traces of derivatives give terms which are zero by the equations of motion. Following 
Polyakov [jT^ we refer to such objects as the 'letters' of the N = A 'alphabet'. We will use 
A to denote the collection of letters. 

We define d{A) of a letter A to be its engineering dimension, so d{(^^) = 1, (i(A") = |, 
d{F^i,) = 2, and each covariant derivative adds one to the dimension. The enumeration of 
letters weighted by dimension gives rise to the following elementary partition function for 
the A/" = 4 alphabet [|T5|,|T7 



zix) = Y: = (2.4) 

This function has a power series which converges for all temperatures < x < 1 and can 
be understood as follows: 

z(x) = ^6^ + 16x^+ 30x^ + 48x^ + 70x^ +■■■. (2.5) 

6{<pi) 16(A-) 24(D^<^-f)+6(F^,) 48(D^A-) 54(D^D,0-r) + 16(D^F,p) 
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There are only 48 -D^A" instead of 64 because 16 components are set to zero by the 
Dirac equation. Similarly there are only 54 components of D^D,j(f)^ because D'^cj)^ is zero. 
Finally, 4 of the 24 components of D^F^p are zero by the equations of motion and another 
4 are equal to zero by the Bianchi identity. 

In what follows we will frequently need to know various partition functions with the 
fermion number operator (—1)''^ inserted. These factors can be easily dealt with by making 
use of the fact that bosonic and fermionic operators respectively have integer or half-integer 
dimensions at tree level. To this end we introduce the quantity 



uj = e^^^ (2.6) 

which is equal to +1 if it is raised to an integer power and —1 if it is raised to a half-integer 
power. To see uj in action consider the formula 

z{i^x) = e^-^^^^a:'^^^) - 2x(3 + v^) ^ _ ^g^3/2 ^ . . _ ^ (2.7) 

2.3. Single-trace operators 

The next step is to string individual letters of the alphabet A together to form single- 
trace operators. The quantity we would like to calculate is the single-trace partition 
function 

Z{x) = x^"(^). (2.8) 

single— trace O 

(Note here the notational distinction between the single-trace partition function Z{x) 
and the full partition function in ( |2.3| ), denoted Z{x).) The only constraint on O = 
Tr[AiA2 ■ ■ ■ A/,] is the overall cyclic invariance of the trace. The enumeration of such 
single-trace operators is identical to the combinatorial problem of counting the number of 
distinct necklaces composed of a collection of different types of beads. The solution to this 
problem, which we review in the next section, may be expressed in terms of z{x) as 

Z{x) = -z{x) - Y ^ - z(^^"+ix")], (2.9) 

n=l ^ 

where (f){n) is Euler's totient function which counts the number of integers less than n 
which are relatively prime to n. The first term —z(x) in ( |2.9| ) is present simply to subtract 
away traces of a single letter Tr[74], since these vanish automatically in the SU{N) theory. 
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Plugging ( p.4| ) into (2^) gives an expansion which goes as foUows: 



Z{x) = 21x^ + 96x2^ 

21(Tr[(t>i(t>J]) 96(Tr[<^-f A"]) 

+ 376x^ + 

76(Tr[<^-f<^J<^-!f]) + 144(Tr[<?!>-rD^0-^])+36(Tr[0-fF^,]) + 12O(Tr[A-A'']) 



(2.10) 



Two comments are in order. The first is that the result ( p.9| ) is only valid in the N = oo 
limit of the the SU{N) gauge theory, since we allow arbitrarily high powers of the individual 
letters. If N were finite, then trace identities would allow a single trace of more than 
letters to be reexpressed in terms of higher-trace operators, indicating that the basis of 
operators we are using would be overcomplete. 

The second observation is that the power series expansion of ( p.9| ) converges only for 
< a; < xyl, where 

= 7 - 4^3 ^ 0.072. (2.11) 

The divergence of the partition function at this value, which corresponds to the tempera- 
ture 

Th(A = 0) = --^ = l/ln(7 + 4V3) ^ 0.380, (2.12) 
mxH 

(measured in units of R~^) has been argued to be the gauge theory dual of the Hagedorn 
transition in string theory [|T5| , p!7[| . Hagedorn-like behavior in other free large systems 



was also observed in [|0 



2.4- From single-trace to multi-trace operators 

Having determined the partition function Z(x) for single-trace operators, it is an easy 
combinatoric problem to calculate from this the partition function for an arbitrary number 
of traces, since the only constraint is that traces should be treated as indistinguishable 
bosons or fermions. The result is 



Z{x) = exp 



^ n 

,n=l 



(2.13) 



The partition function Zk{x) for /c-trace operators can be extracted by inserting y'^ into 
the sum and then reading off the coefficient of in the expansion of the exponential. 
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From the results reviewed above we can see that the complete partition function of 
A/" = 4 SU{N) gauge theory at infinite and zero Yang-Mills coupling is given by the 
formula 



exp 



n 



1 



^2.14) 



The exponential term in (|2.14| ) is (one over) the partition function for gauge group U{1) 



and the infinite product is the partition function for gauge group U{N), so ( |2.14| ) expresses 
the expected fact that 



Z 



SU(N) 



U{N) 



Z, 



c/(i) 



(2.15) 



In fact, it is interesting to note that although this analysis has been done at infinite A^, 
the result ( p.l4|) remains correct to all orders in 1/A^ (though certainly not at finite A^). 
This is true because the tree-level dilatation operator Dq obviously does not receive any 
1/A^ corrections, and trace relations are non-perturbative in 1/A^. 



2.5. Turning on the Yang-Mills coupling 

In (|2.14| ) we have reviewed the complete partition function for SU{N) A/" = 4 SYM 
theory on M x S''^ at infinite and zero Yang-Mills coupling. The goal of this paper is to 
calculate the first-order correction to this result when we turn on the Yang-Mills coupling 
Qym- The dilatation operator D can be expanded in powers of the 't Hooft parameter 

^ = 9ym^ as 

D = Do + -^D2 + ---, (2.16) 
so to first order in A the partition function ( |2.tj| ) is given by 

Z(x, A) = Tr[x^o+i^^^+-] = Tr[x^o] + Al^Tr[x^"L>2] + ■ ■ ■ ■ (2.17) 
Therefore we need to calculate the trace of the one-loop dilatation operator, 

^(i)(x) = Tr[x-^°D2]. (2.18) 

Our calculation will proceed by first calculating the one-loop partition function in the 
single-trace sector and then assembling together the result for multi-trace operators as in 
the previous subsection. 
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2.6. Operator mixing and 1/N 

The tree-level dilatation operator Dq is diagonal in the trace basis, since the engineer- 
ing dimension of a /c-trace operator is obviously just the sum of the engineering dimensions 
of the k individual operators. At one loop this is no longer true. Instead we have a formula 
of the form 

k 

AT' ' A^' 



k 

D2\k) = J2 T^[Oi] ■ ■ ■ {D2Tr[0,]) ■ ■ -TrfOfc] + hk - 1) + hk + 1) (2.19) 



i=l 



for a generic /c-trace operator = Tr[Oi] • ■ ■ Tr[Ofc]. At first glance one might be tempted 
to disregard the 0{1/N) terms since we are working at infinite A^. However it is well- 
known that there are classes of /c-trace operators whose one-loop anomalous dimensions 



receive non-zero contributions from mixing with k ± 1-trace operators pT| , |2^ . This occurs 
generically for BMN operators, which consist of L ~ ^/N letters. For such operators the 
1/A^ suppression is overwhelmed by the growth of the number of non-planar diagrams ||23|| . 

Therefore we will not use A^ = cxo as a justification to omit the last two terms in 
( p.l9|) . Fortunately we have an even better excuse, which is simply that these terms are 
non-diagonal in the trace basis and therefore do not contribute to the quantity ( |2.18| ) 
that we are computing. The situation becomes more complicated starting at two loops, 
where the correction to the partition function involves Tr[x^°I^4] and Tr[x^° 02], both of 
which have some diagonal terms of order l/N"^ which can't necessarily be dropped. (In 



some sense this smells like an 'order of limits' problem, like the one discussed in ||2^). Of 



course, if one focuses on calculating the planar partition function (as opposed to the large 
A^ partition function), then none of these issues arise and one can ignore the 1/A^ terms in 
( p.l9|) from the beginning. Subtleties in the 1/A^ expansion have been noted in [p5| . 

Note that we have implicitly been using a scalar product on the space of operators 
which is diagonal in the trace basis, 

{k\l) oc 5ki. (2.20) 

The Af = 4 SYM theory provides, via the state-operator correspondence, a natural inner 
product Ski = {k\l) on the space of operators which is diagonal in the trace basis at infinite 
A^ but receives non-diagonal, trace-mixing corrections beginning at 0{1/N). This operator 
mixing does not concern us here since the trace 

TT[x''°D2]=Y,{k\x'''D2\l)S-^ (2.21) 

k,l 
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is completely independent of the scalar product S. Therefore we are free to choose the 
most convenient inner product Ski = ^ki- A non-trivial inner product S appeared in studies 
of corrections to the BMN correspondence, such as [p6| , |27[1 where matrix elements of 
the Hamiltonian were compared between gauge theory and string theory. The utility of 
ignoring the gauge theory inner product for the purpose of calculating basis-independent 
quantities was emphasized in in the context of calculating eigenvalues of D2 in the 
BMN sector. 



3. Polya Necklaces 

In this section we develop the machinery which will reduce the calculation of the 
one-loop partition function 

Z^^\x) =Tv[x^°D2] (3.1) 

in the single-trace sector to the calculation of some elementary traces involving O = D2 
acting on two neighboring letters. Roughly speaking, we consider here the problem of 
'tracing out' all but two letters of any gauge-invariant operator, and express the result in 
terms of traces of O over the remaining two letters. This is accomplished by translating 
the calculation into the language of spin chains, and then using a generalization of Polya's 
counting theory. 

3.1. Free necklaces 

A necklace of length L is a collection of L objects [Ai ■ ■ ■ Al) chosen from some 
fixed set A E A of 'beads', such that two necklaces are identified if they differ from each 
other by a cyclic rotation. It is useful to introduce a counting function d on the beads, 
and define d to act additively on the beads of a necklace, 

L 

d{A^---AL) = Y,diA,). (3.2) 

1=1 

The analysis of this section will be general, but of course for the desired application to 
A/" = 4 SYM theory we remember in the back of our mind that we will take A to be the 
A/" = 4 alphabet, and d{A) will be the engineering dimension of the letter A. 

A central result of Polya's counting theory is that the necklace partition function 

Z(x) = 5^x'^W, (3.3) 

N 
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where we sum over necklaces of arbitrary length L > 1, is given by 

Z(x) = -f;^ln[l-^(x")], (3.4) 
n 

n=l 

where z{x) is the generating function for the beads, 

^(x)=^x^(^). (3.5) 
AeA 

The formula ( p.4| ) is valid when the beads are all bosons. The generalization to include 
fermions is straightforward and will be presented below. 

3.2. Necklaces with a pendant 

Instead of reviewing the elementary derivation of the result ( p.4|) , we will consider a 
useful generalization from which we will recover ( |3.4]) as a special case. To describe the 
generalization that we are interested in, it is useful to think of a necklace of length L 
as a spin chain of length L, where on each site of the spin chain the spin vector takes 
values in the set A. Of course we have the constraint that only cyclically invariant spin 
configurations correspond to necklaces. Therefore, we can recast the calculation of ( |3.3| ) 
into spin chain language by writing the partition function as 

oo 

N L=l 

where Tri, denotes the trace over spin chains of length L and V denotes the projection 
operator onto the subspace of cyclically invariant spin configurations. The projector can 
be written explicitly as 

V = y{l + T + T^ + ---+T^-^), (3.7) 

where T is the translation operator which sends site i on the chain to site i + 1 and satisfies 
= 1. 

The generalization of TvLi'Px'^] that we would like to consider is to the case 

TrL[Vx''di (3.8) 

where O is any homogeneous operator which commutes with d and which acts only on two 
neighboring sites at a time, so that it may be decomposed into the form 

L 

^ = J20^,^+l, Oi,i+i = 1i®---®1,_i®0®1,+2®---®1l. (3.9) 

2=1 
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Spin chain Hamiltonians are of course prime examples of such operators, and eventually 
we will apply the present machinery to the case O = D2, but our analysis will continue 
to be as general as possible. Since the operator O only connects two neighboring sites at 
a time, but can slide all the way around the necklace, we can think of O as a 'pendant' 
hanging from two adjacent beads on the necklace. 

Given any such operator O, an obvious quantity of interest is the expectation value 

{0{x))=TTy,^A[x''0]= Yl a;"^^^^+'^^^^^i^2|0|AiA2). (3.10) 

We will see that knowing (0(x)) is almost, but not quite enough information to allow for 
the calculation of ( |3.8| ). The other quantity that we will need to know is the 'permuted 
trace' 

{PO{w,y)) = TrAxA[Pw'''y'''0]^ u;'^(^^)/(^^)(AiA2|0|A2Ai), (3.11) 

Ai,A2eA 

where P is the permutation operator on ^ x ^ and the two sites are counted with different 
variables w and y. Let us now see how to reduce the calculation of ( p.8| ) to the calculation 
of these two quantities. 

Since V projects onto cyclically invariant states anyway, the sum over i in (|3.9| ) is 
actually redundant. We can affix the pendant to sites (i,i + 1) = (1,2), and the sum in 
( P^) just gives a factor of L which cancels the 1/L in ( p.7| ) giving 

L-l 

Tr^iPx^O] = TrL[T''x'^di2]. (3.12) 

k=0 

This trace may be expressed as 

L-l 



TTL[Vx''d] = Y J2 {Ai---AL\T''x''di2\Ai---AL) 

k=0 Ai,...,Al 

k=0 Ai,...,Al i=3 



(3.13) 



Upon contemplating the formula (|3.13|) , it is clear that after we sum over A^, . . . , A^, 



only two possible structures can emerge, 

(z) (A1A2IOIA1A2), or (n) {AiA2\0\A2Ai), (3.14) 
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depending on whether k and L are such that the Kronecker delta functions in (|3.13 ) end 



up connecting site 1 + /c to site 1 or to site 2. For example, k = clearly gives the former 
structure while k = 1 clearly gives the latter. But for general k and L, what is the criterion 
which tells us which of the two possibilities ( |3.14|) we get? 

Let us start with site 1 + k and follow it around the necklace using n of the Kronecker 
delta functions, 

1 + A; ^ 1 + 2A; ^ 1 + 3A; ^ >l + {n + l)k. (3.15) 

Since there are only L—2 delta functions in total, we have < n < L—2. Now let m = (L, k) 
be the greatest common divisor of L and k. If m > 1, then by step n = L/m — 1 at the 
latest, site 1 + k will have been connected to site 

l + A;^l + (n + l)A; = l + (L/m - 1 + 1)A; = 1 + {k/m)L = 1 mod L. (3.16) 

On the other hand, if k and L are relatively prime (m = 1), then there is no n < L such 
that 1 + (n + l)k = 1 mod L, so 1 + k cannot get connected to site 1 and hence must be 
connected to site 2. We conclude that we get the trace structure of type (m) if and only 
if k and L are relatively prime, and we get structure (z) otherwise. Let us therefore study 
these cases separately. 

3.3. Case (i): k and L have a common factor 

As a warm-up exercise let us consider the case L = 15, /c = 6, so that m = (15, 6) = 3. 
Then as we sum over sites 3, 4, ■ ■ ■ , L, the delta functions in ( |3.13|) sew together the beads 
of the necklace as follows: 

7 ^ 13 ^ 4 ^ 10 ^ 1, 

8^ 14^5^ 11 ^2, (3.17) 

9 ^ 15 ^ 6 ^ 12 ^ 3 ^ 9. 
Each line in this formula can be though of as a 'strand' of the necklace. So the necklace 
with L = 15 and /c = 6 is composed of m = (15, 6) = 3 strands, and the Kronecker delta 
functions in ( |3.13D force all of the beads on any given strand to be the same. 

The strands starting with 1+k and 2+k end up respectively at sites 1 and 2, confirming 
our earlier analysis. The third line in ( |3.17| ) denotes the following contribution to ( ^.13| ): 

A3,AQ,Ai5,Ae,Aj_2 

A3 

(3.18) 
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The first two strands in ( p.l7|) involve the sites 1 and 2 where the pendant is attached, 
and it is not hard to see that they end up contributing a factor of 



J2 x^^(^i)+^'^(^2)(AiA2|0|AiA2) = (0(x^)), 

Ai,A2 



(3.19) 



using the definition ( p.lO| ) . Combining these results, we find that the total contribution to 
the sum ( p.l3| ) for the case L = 15 and /c = 6 is 



z{x^) (0(x^)). 



(3.20) 



The generalization of this analysis is straightforward. A necklace with general L and 
k will have m = (L, k) different strands, each with L/m beads. Two of those strands (like 
the first two in ( p.l7|) ) will involve the sites 1 and 2 and give rise to 



(0(x-^/"^)). 



(3.21) 



The remaining m — 2 strands (note that we are assuming here that m > 2), like the last 
line in ( |3.17| ), give a factor of 



Combining ( p.21|) and (|3.22|) , we conclude that the total contribution to (|3.13|) from all k 
such that m = (/c, L) > 1 is 



L-l 



k = 
m=(k,L)>l 



r -1 m — 2 

J2 k^^/"") (o(x-^/'^)). 



(3.23) 



3.4- Case (ii): k and L are relatively prime 



As a warm-up exercise let us consider the case L = 14, k = 5. The Kronecker delta 
functions in (|3.13|) now sew together the beads 



6 ^ 11 



7 ^ 12 



8 ^ 13 



9 ^ 14 



10 ^ 1, 



(3.24) 



confirming our general analysis that site 1 + k gets connected to 2, and site 2 + k gets 
connected to 1, giving trace structure [ii). The first line indicates a contribution of 



X 



d(A6)+d(Aii)+d(A2) ^ ^3d(A2) 



(3.25) 



Ae.Aii 
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while the second Une similarly denotes a contribution of x^^^*-^^-*. Putting everything 
together, we find that for L = 14 and k = 5, ( |3.13| ) reduces to 



J2 x^^'^^^'^+^'^^^'\AiA2\0\A2Ai) = {POix^\x^)), (3.26) 



Ai,A2 



using the definition ( |3.11| ). 

The generalization to arbitrary k and L is straightforward. After n — 1 steps in the 



first line of (|3.24|) , site 1 + k will connect to site 1 + nk. Therefore, the length of the 
first strand is the smallest non-negative n such that 1 + nk = 2 mod L, or equivalently 
nk = 1 mod L. The total contribution to (|3.13| ) from all k which are relatively prime to L 
is therefore 

L-1 

J2 (PO(x-^-"(^''=), x"^^''^))), n(L, k) = mm{n > : nk = 1 mod L}. (3.27) 

fc=0 

(fc,i)=l 

Solving nk = 1 mod L for fixed k and L is equivalent to finding n, m such that nk—Lm = 1. 
Given that k is relatively prime to L, it is clear that a solution exists only if n and L are also 
relatively prime. Moreover, it is clear that n < L (otherwise subtracting L from n would 
give a smaller solution). Finally, the set of {k : {k, L) = 1} is in one-to-one correspondence 
with the set of n(L, k), simply because the condition nk = 1 mod L is symmetric in n and 
k. Therefore, although n(L, k) is not generically equal to /c, the sum in (|3.27| ) is equivalent 
toS 

L-1 

J] (PO(x^-^x'=)). (3.28) 

fc=0 
(k,L) = l 

3.5. Summary and main result 

We now combine the contributions (|3.23| ) and (|3.28| ), writing the result as 

i{k,L)-2 



fc=0 

L-1 

+ [(PO(x^-^x'=))-^(x^)-l(0(x^))] 



(3.29) 



fe=0 
(k,L) = l 



^ We are grateful to Jan Plefka for pointing out a flaw in our original proof of this formula. 
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where in the first term we omitted the constraint m > 1 from ( p.23|) at the expense of 
subtracting off the extra terms on the second line. Now we can trade the sum over k in 
the first term of ( |3.29| ) for a sum over divisors a of L, to write 



L-l 



fc=0 
(k,L) = l 



(3.30) 

At this step let us pause for a moment to explain how to recover the Polya formula 
( ^.4| ) as promised. To this end we need to consider the special case where the operator 
O is proportional to the identity matrix, and specifically O = 1/L. This is the correct 
normalization which gives rise to O = 1 when plugged into the sum over sites in (|3.9[) . For 
O — 1/L we easily find 



(0(x)) = {PO{w, y)) = ^{wy). 



The second term in (|3.30|) therefore drops out, leaving just 



(3.31) 



(3.32) 



The sum over L is performed in the usual way, and we obtain 



TtlIVx"] = _ ^ M in[l - z{x 

Til 



L=l 



n=l 



(3.33) 



which is the desired result (|3.4| ). 

Having confirmed that the formula ( p.30|) reduces to the known answer for the special 
case O = 1, let us now consider operators O which do not depend explicitly on L. In 
particular this implies, through ( ^.9] ), that the eigenvalues of O scale linearly with L. 
Then we can perform the sum over L > 1 to arrive at 



^Tri[Px^O] 



L=2 



z{x) 



+ ^(/>(n 



(0(a;'^)) 



1 



n=l 



z{x'^) 1 — z{x'^) 



(k,L) = l 



(po(x^-^x'=)) - 



z{x^) 



(3.34) 



The first term is present to subtract off the part of the second term which would correspond 
to L = 1, which we omit since it can't support a pendant (and moreover is irrelevant in 
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the SU{N) gauge theory). A final step is to simphfy ( p.34|) by changing the summation 
variable from n to L and combining everything into the main result 



L=2 L=lk:{k,L) = l 



1 - z{uj^^^x^) 



(3.35) 

It is a straightforward exercise to generalize the analysis of the previous subsections to 
allow for fermionic beads, and we have included here the appropriate factors of uj which 
keep track of the minus signs appearing when such beads are permuted. (The permutation 
operator P is understood to be graded, i.e. P\AiA2) = (— 1)^^^^ IA2A1).) 



The first term in ( p.35 ) for L = 1 is precisely what one would obtain by making the 
crude estimate that the only effect of the projection V onto cyclically invariant states is to 
insert a factor of 1/L. The rest of ( |3.35| ) is the detailed correction to this approximation. 

In all the cases relevant to A/" = 4 SYM theory that we study below, we will see that 
the second term in ( |3.35|) is a very small correction in the sense that its contribution to 
the coefficient of x"^ is negligible for large n. In particular, we will find that {0{x)) and 
{PO{w, y)) converge for all temperatures so that the large temperature behavior of ( |3.35| ) 
is dominated by the pole 1/(1 — z{x)) in the first term. 



4. Examples 

We can gain some insight into the formula (|3.35| ) by applying it to some subsectors 
of the gauge theory. The implication of the results presented in this section for the ther- 
modynamics of A/" = 4 SYM theory is unclear since there is no sense in which the sectors 
decouple from each other at finite temperature (we do not consider here the addition of 
chemical potentials for various charges) . However, we believe this section is a useful warm- 
up exercise for the more complicated analysis which follows. Moreover, the results given 
here for the traces of the SU{2), 5*0(6) and SL(2) spin chain Hamiltonians may be of 
interest from the point of view of integrability in those sectors. An additional subsector 
which is of independent interest is the SU {2\4) subsector, which at one loop is isomor- 



phic to the 't Hooft large limit of the plane- wave matrix model [10]. This subsector is 



considered in [E9 . 
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4.1. The SU{2) subsector 

This subsector consists of all operators of the form 



Tt[XXZZZXZZ-XXZ], (4.1) 

where X and Z are two holomorphic scalar fields. The 'alphabet' for this sector is ^ = 
{X,Z}, the dimension formula is d{A) = 1, and the elementary partition function is 
z{x) = 2x. In the natural basis {XX, XZ, ZX, ZZ} for Ax A, the matrix elements of the 
permutation operator P and the one-loop Hamiltonian D2 are simply 



/I 0\ 

10 

10 

Vo 1/ 



D2 = l{l-P). (4.2) 



We immediately find 

{D^ix)) = x\ {PD2{w, y)) = -wy. (4.3) 



Plugging these into the main formula ( p.35| ) gives, after some simplification, the following 



formula for the trace of the Hamiltonian in the SU{2) sector: 



Trix^'-D^] =x-J2 <Piri)x- ( — — ^ ) . (4.4) 



n=l 



As a check, we used a computer to calculate the trace of the SU{2) spin chain Hamiltonian 
for all chains of length L < 26 and successfully matched the coefficients in the expansion 
of (|4.4|) up to order x'^^. 

4.2. The SO{6) subsector 

This subsector consists of all operators built only out of scalar fields with no deriva- 
tives, 

Tr[cl)'^---cl)'^], /, = {1,...,6}. (4.5) 

The alphabet is ^ = {^^, . . . , (/)^}, the dimension formula is d{A) = 1, and the elementary 
partition function is z{x) = 6x. In the natural basis I/1/2) = (f)^^(f)^'^ for Ax A, the matrix 
elements of D2 and P are 

{Ill2\D2\JlJ2) = -^Si^ijj^j^ + -Sl^Jj^Si^j^ - -Si^jji^j^, 
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A simple calculation yields 



{D2{x)) 



33 



-X 



27 

{PD2{w,y)) = -^wy, 



(4.7) 



which leads to the result 



27 3 ^ 

— X > 

2 2 ^ 



71 X 



9 - 65x^ 
1 - 



(4. 



for the trace of the 5*0(6) Hamiltonian. As a check, we used a computer to calculate the 
trace of the SO{Q) spin chain Hamiltonian for all chains of length L < 11 and successfully 
matched the coefficients in the expansion of ( ^?8|) up to order x^^. 

4-3. The SL{2) subsector 

This subsector consists of all operators of the form 



Tr[D''Z---D'^Zl in E {0,1,...}, 



(4.9) 



where Z is a single holomorphic scalar field and D is a single holomorphic covariant deriva- 
tive. The alphabet is ^ = {Z,DZ,D^Z, ■ ■ ■}, the dimension formula is d{D^Z) = i + 1, 
and the elementary partition function is 



4x) = 



X 



d{A) 



AeA 



X 



i=0 



X 



(4.10) 



In the basis \i1i2) = D^^Z D^^Z for Ax A the matrix elements of D2 are ||30|| 



{iii2\D2\jij2) = 2'^M+*2,ji+j2 



where h{j) are the harmonic numbers 



(4.11) 



Mj) = E^, Mo) = o. 

n=l 

The matrix elements of P are obviously 



(4.12) 



{ili2\P\jlj2) = 5^^JJ^^j^. 



(4.13) 
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A simple calculation yields 



X 

(1-X)2 



ln(l — x) 



(4.14) 



and 



2 1 — wy 



(1 - wyY 



(4.15) 



After some simplification, we find for the trace of the SL[2) Hamiltonian the result 



X^ 



IY[x^oi^2] = -5^</>(n)^^^ln(l-x-) + 5^^4^ In(l-x^). (4.16) 



n=l 



L=l 



fe = i 
(fc,i) = i 



As a check, we used a computer to calculate the trace of the SL{2) spin chain Hamiltonian 
for all chains with total dimension Dq < 19 and successfully matched the coefficients in 



the expansion of (|4.16|) up to order 



X 
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5. Traces of the PSL{4\4) Spin Chain Hamiltonian 

We now turn to our next step, which is to apply the result ( p.35| ) to the calculation of 
the one-loop correction to the partition function of A/" = 4 SYM theory on M x S""^ in the 
single-trace sector: 

Z^^\x) =Tr[x^''D2]. (5.1) 
To this end, we compute in this section the quantities 

{D2ix)) = T:^x^[x^°I^2], {PD2{w,y)) = TrAxA[Pw''^'(^^y'''^^'W2] (5.2) 

needed to invoke ( |3.35| ) for the full P5'L(4|4) spin chain Hamiltonian D2. 

The calculation of {D2{x)) is greatly facilitated by making use of the P5'L(4|4) symme- 
try of A/" = 4 SYM theory. Since the dilatation operator D2 commutes with this symmetry, 
the action of D2 on an arbitrary state can be decomposed into its action on irreducible 
representations of P5'L(4|4). We therefore begin with a discussion of the relevant repre- 
sentations. Unfortunately, the operator w^'''-^''y^°^'^'> does not commute with the two-letter 
PSL{4\4) Casimir, so the calculation of {PD2{w,y)) will prove more difficult. 
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5.1. Digraphs in the Af = 4 language 

The elementary fields and their covariant derivatives which make up the alphabet A 
of A/" = 4 SYM theory constitute the so-called 'singleton' representation of PSL{4\4). The 
superconformal primary state is the scalar field (f)^ , with quantum numbers 

[0,l,0](o,o) (5.3) 

under SL{4) flavor rotations and the SL{2) x SL{2) Lorentz algebra. The singleton rep- 
resentation is frequently denoted Vf, although we shall continue to refer to it as A for 
consistency. 

Since the one-loop dilatation operator D2 only acts on two letters at a time and 
commutes with the two-letter Casimir of PS'L(4|4), it is sufficient to consider the 

decomposition of the product of two copies of the singleton representation into irreducible 
representations of PSL{4\4). The decomposition is 

00 

AxA = Y,^j^ (5-4) 

where Vj is the module whose superconformal primary is an eigenstate of the P5'L(4|4) 
Casimir with eigenvalue j{j + 1) and quantum numbers 

[0,2,0](o,o), [l,0,l](o,o), and [0, 0, 0](^_i^^_i) for j > 2. (5.5) 



In the notation of ||3T[|, we have 



A = bI'} 



qo,i,o](o,o)' 



1 1 

'[0,2,0] (0,0)' 



(5.6) 



[1,0,1](0,0)' 

andV,=Cf,;^„_„,(.^._^^.^._,)forj>2. 



In linguistics, a group of two successive letters whose phonetic value is a single sound, such 
as ng in Yang or th in theory, is called a 'digraph', so we can think of the Vj as the digraphs 
of A/" = 4 Yang-Mills theory. 

It is straightforward to count the primary states in Vj, weighted in the usual way by 
. For j = and j > 2 the results can be read off respectively from (6.13) and (5.45) of 
31j] , or from tables 7 and 8 of [ffSl . We did not immediately find the primary content of Vi 
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in the literature, but the derivation thereof is straightforward and the result is presented 
in appendix A. The results for all Vj can be summarized in the expressions 



X 



Vo{x)=4x'{l + ,/^)\5 

(5.7) 

Vj{x) =x\l + V^yU - 1 + (i + 2)^){j - 1 + 5v^ - (j + 2)x), j > 1. 

Setting X = 1 counts the total number of primaries in Vj, which is 2^(2j + 1) for any j > 0. 
Note that for j = 0, 1 some powers of x in Vj{x) have negative coefficients. This may 



be thought of as a bookkeeping device (explained in [|31[]) which allows for easily keeping 
track of fields which are eliminated by the requirement of imposing equations of motion 
or conservation laws. One consequence of this is that the full partition function for the 
module Vj (including descendants) may be calculated naively, with derivatives treated as 
if they acted freely, without worrying about equations of motion or conservation laws. The 
partition function for the module Vj is therefore simply 

lYv, = Tr^x^[P,x^"] = (5-8) 

where we have defined Pj to be the projection operator Pj : ^ x ^ ^ Vj. 
The Pj form a complete set of orthogonal projection operators, so 

oo 

5^P, = 1. (5.9) 



This implies the identity 



Y.T^AxAiPjx''''] = IV^x^[x^°] = z{x)^ (5.10) 



which is indeed satisfied by (|5.7|) and (|5.8D . 

5.2. Simple trace {D2{x)) 

Here we calculate the expectation value ( |3.10| ) for the one-loop dilatation operator 
D2. The calculation only takes one line since we have all of the machinery in place. In 



it was shown that the eigenvalue of D2 in the module Vj is simply the harmonic number 

Mj) = 5Z-, /i(0) = 0, (5.11) 



n=l 

21 



and therefore that D2 may be written as 



00 



D2 = Y.h{j)P,. 



(5.12) 



j=o 



Prom (|5.8| ) and ( |5.12| ) we immediately have 



00 




(5.13) 



Plugging in ( p.7|) and performing the sum over j gives 



(1 + y/^)^ 
(1 - v^)6 



[-(1 - 4v/x + xf ln(l - x) - x(l - 8v^ + 2x)] . 



(5.14) 



5.3. Permuted trace {PD2{w,y)) 

As mentioned above, the calculation of {PD2{w, y)) is complicated by the fact the two- 
letter PSL{4:\4) Casimir operator J'^12) commutes only with the sum -Do(i) + -^0(2) 5 but 
not with -Do(i) ^ind -Do(2) separately. In particular, it does not commute with w^°'-^'^y^°'-'^\ 
so the beautiful decomposition into the modules Vj that we employed in the previous 
subsection is of no use here. This apparent breaking of P5'L(4|4) is merely an artifact of 
our choice to simplify the calculation of Tyl\Px^'-'' 02] by tracing out L — 2 sites on the 
chain and expressing ( |3.35| ) in terms of the two remaining sites. 

A manifestly P5'L(4|4)-invariant calculation ofTrL[Vx^° D2] would proceed as follows. 
First, one would need to know the collection Cl of PSL{4\4) modules which appear in 
higher powers of the singleton representation A, 



(where Tr denotes the projection onto singlets of the cyclic group I^l)- Then PSL{4\4) 
invariance guarantees that in each resulting module Vj the dilatation operator D2 is pro- 
portional to the identity operator, with some calculable eigenvalue hj. The desired trace 
would then be given by 




(5.15) 



oc 




(5.16) 



L=2XeCL 



where Px is the projection operator from Tr(^^) onto Vx- 
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The decomposition of ^ x ^ x ^ appears in (see also 0), where it was used 
to determine the one-loop anomalous dimensions of some operators consisting of three 
elementary fields. However, it seems quite challenging to implement the procedure outlined 
in the previous paragraph for arbitrary L, although of course it would be very interesting 
to do so. 

Instead, we will proceed by using the matrix elements of the operator L>2, written 



down in in a GL{4\4) oscillator basis, and then calculating the desired trace 'by hand'. 
This is quite a lengthy calculation, so we begin by presenting the result. The interested 
reader can find more details below. We find: 



{PD2{w,y)) 



1 



(1 - v^)3 (1 - ^)3 (1 + ^)3(^ + ^) 



wy{l - w){l -y)pi{w,y) 



1 , 

+ - In 
2 



{l-w){l-y) 
(1 - wyY 



(1 - ./wy){^/w + ^/yfp2{w,y) 



1 iVw- ^) 

2 (v^ + v^) 



In 



1 — w 



(1 + ^/wyfpsiw.y) 



in terms of the three polynomials 



(5.17) 



pi = 4 - 16w^/^ + 7w- 16y^/^ + 22w^/^y^/^ - 16wy^/^ 

+ w^/^y^/^ + 7y- IQw^^^y + Qwy + w^^^y^/^, 



(5.18) 



P2 



1 - 4w^/^ + w- 4yi/2 + 20w^/^y^/^ - 20wy^/^ + 4w^/^y^/^ +y- 20w^/^y 

+ 42wy - 2Qw^'^y + w'^y + 4w^^^y^''^ - 20wy^/^ + 20w^/'^y^/^ (5.19) 

A 2 3/2 I 2 A 3/2 2 I 2 2 

— 4w y ' + wy — 4w ' y + w y% 



and 



P3 = w - 4w^/'^ + w'^ + 4w^/^y^/^ - 20wy^/^ + 20w^/^y^/^ - 4w^y^/^ + y 

- 20w^/^y + 42wy - 20w^/^y + w^y - 4y^/^ + 20w^/^y^/^ (5.20) 

- 20wy^/'^ + 4m;3/2^3/2 + ^2 _ 4^1/2^2 ^ ^^2_ 

Note the interesting identity 



P3{w, y) = y^p2(l/y, w) = w'^p2{y, l/w). 



(5.21) 
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It is useful to subject this complicated result to a simple check. When we set w = y = 
X, then the calculation of {PD2{x,x)) can be done using the group theoretic techniques of 
the previous subsection. In particular, we have 

{PD2{x,x)) = TrAxA[Px'''D2] 



^Trv,[Pa:^°I)2] 

j=0 



(5.22) 



Vj{x) 



j=0 



(i-xY 



where we used the fact that the permutation operator P acts as ( — 1)-' in Vj [ p3| , p3| . If we 
now substitute the expressions for Vj{x) from ( |5.7| ) and perform the sum over j, we obtain 

+ ^ — \xil - 7x^2 +X + x'/' - 6x^ + 2x^/2) 

(5.23) 



{PD2{x,x)) 



(l-v^)4(H-^)'L ^ ' ^ 

- (1 - 7x^/^ + 15x - 25x^/^ + 25x2 - 15x^/2 + - x^/^) ln(l + x) 

Encouragingly, this expression agrees precisely with that obtained by setting w = y = x 
in (^. 

Now we begin in earnest the calculation of ( |5.17| ). The first step is to use the matrix 



elements of D2 in a GL(4|4) oscillator basis, as presented in ||3y], to write down a sum 
which gives (|5.17[) : 

00 4 
{PD2{w,y))= Yl E (-1) 



si,S2,Pi,P2,k=0 Fi,F2,j=0 



4-j\f4-j 



jj \F2-jJ k\ik + l)\ 



X c(ni + n2,ni-k- j, n2 - k - j)^yi+-i/2+pi/2^i+s2/2+p2/2 



X - k, -k, -si, -S2; 2, 1 - /c + 1 - /c + P2; 1) 

2 

X 

i=l 



(5.24) 



Jj5(l- is, + - iF,)(l + s,)(l+Pi 

i=l 

where J-" is the regularized hypergeometric function, 

Hi = Si + Pi + Fi, 

and the coefficients c are the matrix elements of D2 given in 

c(n,ni2,n2i) - -(-1) r(l + W) 



(5.25) 



(5.26) 



c(n, 0,0) = i/i(in). 
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The detailed derivation of ( ^.24 ), which is not entirely straightforward, is presented in 
appendix B. The next several steps of the calculation will be shown schematically. The 
skeptical reader is free to verify that the power series expansion of ( p.24| ) agrees with that 
of (|5.17|) to any desired order. 

After several manipulations similar to the ones in appendix B, equation ( p.24|) can be 
cast into the form 



{PD,{w, y)) = Y^i-inn + l){n + 3)w'+''/\ 



i/2 



n=0 
X 



(1 - vyy 



(5.27) 



where Pniv) is a polynomial in y^y whose highest term is 0{y"'~^^). We did not obtain 
an explicit formula for Pn{y), though presumably it could be reverse engineered from the 
final answer ( |5.17| ). Instead, we use a trick by investigating the quantity 



Q[w,y) = {PD2[w,y)) ^ 

L(l + v^) (1 + vw) 

Now we use ( |5.27 ) and the power series expansion 

• -1 / — \ 2 oo 
1 - y/w 



(5.28) 



1 + ^yw 



m=l 



i/2 



dm = Sm,o + 4m(-l)' 



(5.29) 



to write 



Qiw,y)= J2 dm{-inn + l){n + 3)w'+^/'+^/^{y-^/'P^iy) + W^iy) 



m,n=0 



where 



WM=y- 



i/2 



(l + (_l)-+lj/l+-/2)2 



Hl-y) 



(5.30) 



(5.31) 



The trick is now to break the power series expansion oiQ{w,y) into the upper diagonal 
terms (where the power of y is greater than the power of w), the diagonal terms, and the 
lower diagonal terms. Since Q{w, y) is a symmetric function, the lower diagonal terms will 
be known once the upper diagonal terms are known. Furthermore, the diagonal terms can 
be extracted from (|5.23| ), so all we have left to calculate are the upper diagonal terms. 
But since Pn is a polynomial whose highest-order term is y^~^^, we see from (|5.30| ) that it 
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never contributes to the upper diagonal. Therefore, for purposes of computing the upper 
diagonal we are free to omit the Pniv) term in (|5.3UD, and make the replacement 



Wn{y) Wn{y)\yk,k>l+n/2+m/2 



(5.32) 



where the notation means that we write out Wn{y) as a series in y and throw away all 
powers of y less than or equal to 1 + n/2 + n/2. This finally gives a formula for Q{w,y) 
which is amenable to a calculation in Mathematica. In this manner, we obtain after a 
tedious but straightforward calculation the result (|5.17D. 



6. One-Loop J\f = 4 SYM Partition Function 
6.1. Single-trace 

The complete one-loop correction to the partition function of TV^ = 4 SYM theory on 
M X S''^ in the single-trace sector is given by substituting (|5.14|) and ( ^.17|) into (|3.35| ) . We 
will not rewrite the formulas because no significant simplification seems to occur. Instead, 
let us note that the result has the expansion 



^Trix^'W^] = ^ 3x2 ^ 4g^5/2 ^ 384^3 ^ 2064x^/2 + ■ 
47r2 ^ ^ 47r2 



(6.1) 



(Each coefficient receives contributions from both kinds of traces, (-D2) and {PD2).) The 
first term in (|6.1|) encodes the one-loop anomalous dimension of the Konishi operator, 
3/47r2. The second term is 3 x 16, coming from 16 descendants of the Konishi operator. 
The third term is 384 = 6x2-1- 20 x3-1- 104 x 3, which come respectively from the 6 
primary states Tr[(j)^ (j)^ (j)"^], which have anomalous dimension 1/2%"^ according to table 3 
in [Roij the 20 Konishi descendants of the form Tr[(/)''^[(/)"', ^^]], and finally 104 Konishi 



descendants which are traces of two elementary fields. 
6.2. Multi-trace 

Now let us extend the result of the previous subsection to the complete one-loop cor- 
rection to the partition function, including multi-trace operators. As discussed in section 
2.6, the diagonal matrix elements of the one-loop dilatation operator act additively on k- 
trace operators. Therefore, to go from the single-trace partition function to the multi-trace 
partition function we can still use the formula (|2.13[) . Substituting 



Z(x) = zW(x) + ^zW(x) (6.2) 
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and expanding to first order in A, we find that the first order correction to the multi-trace 
partition function is 



(6.3) 



(a factor of 1/n is canceled by Inx^ — nlna;) where we recall that the tree level result 
is written in (|2J^ ). 

Now let us plug the result from ( |3.35| ) into (|6.3| ). The term proportional to (P-D2) 
gives the sum 



n=l L=2k:{k,L) = l a,b=l 

(6.4) 

To see why this equation is true, pick any positive a and 6, and look at the left-hand side to 
see how many times (if any) the term {PD2{u>"'~^^x"' ,u''~^^x^)) appears. This is equivalent 
to asking how many solutions, for given a and b, there are to the equations 



nk = a, nL = a + b 



(6.5) 



for L >2 and k : (k,L) = 1. The answer is that there is always precisely one solution: 
n = (a, 6), L = {a + b)/n and k = a/n. Certainly if n were not the greatest common 
divisor of a and b but some smaller common divisor, then (|6.5| ) would still give solutions 
for k and L, but these would not satisfy the constraint (/c, L) = 1. 

Next we plug the {D2{x)) terms from ( |3.35|) into ( |6.3|) , which gives the sum 



°° °° m (, ,nL+l^nL\\ 



n=l L=l 



k=l 



L\k 



1 - z{u'^+^x'^) ~ ^ 1 - z{u'^+^x^) ' 

(6.6) 



Combining (|6.4|) and (|6.6| ) into (^]^) gives the final result 



Z^^\x) 



A In a; 



k=l ^ ' k,m=l 



(6.7) 



as advertised already in (|1.1|), for the one- loop partition function of = 4 SYM theory 
on M X 5"^, expressed in terms of the free partition function Z^^\ the elementary partition 
function (|2^ , and the traces ( |5.14|) and ( p.l7|) . 



27 



6.3. One-loop Hagedorn temperature 



The partition function Z{x) has a simple pole at the Hagedorn temperature 



Xu — X 



(6.8) 



where c is an irrelevant overall numerical coefficient. To compute the one- loop correction 
Sxu to the Hagedorn temperature, we simply expand 



xn + Sxu — X xyl — X 



1 



Sxu 
xu — X 



+ 



(6.9) 



When we compare this to ( |6.7| ) and recall that xh is such that z{xu) = 1, we note that 
only the k = 1 term in the sum of {D2{x^)) contributes to the double pole at the Hagedorn 
temperaturejl Reading off the residue of this pole, we find 



(5xh = — lim 



X — >xh 



xh - x^ 



Alnx (1^2 (x)) 
47r2 1 — z{x) 



2 AlnxH , 
3^h^^(I^2(xh)). 



(6.10) 



Remarkably, we find from (|5.14|) that {D2{xy{)) =3/4, which gives 



(5xh AlnxH 
xh 87r^ 



and hence 



1 dxu X 



InxH xn 87r^ 



(6.11) 



The one-loop Hagedorn temperature is therefore 



A 



Th = Th(0)(1 + ^ + 



7h(0) 



ln(7 + 4^3) 



(6.12) 



It is encouraging that the sign is positive, consistent with the simple guess that the Hage- 
dorn temperature is a monotonically increasing, smooth function of A from Th = Tq at 
zero coupling to the AdS/CFT prediction that Th ~ A^/'* at strong coupling. 



^ It is possible that the second term in (6.7) develops a pole at the Hagedorn temperature 
after evaluating the sum, in which case it would add a finite correction to the one-loop Hagedorn 
temperature. However, this does not happen in any of the subsectors that we studied, and 
numerical evidence suggests that it does not occur here either. 
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7. Discussion 



In this paper we have presented, in (|1.1| ), the one- loop correction to the partition 
function of SU{N), A/" = 4 SYM theory on M x S"^ at infinite and below the Hagedorn 
temperature. 

Several recent papers including ||T7| , ^^ , |37[] address the thermodynamics and phase 
transition structure of weakly coupled gauge theories with the goal of smoothly connecting 
onto the strong coupling predictions implied by the AdS/CFT correspondence. String 
theory in AdS^ x has both a Hagedorn transition at Th ~ A^/^ and a Hawking-Page 
transition involving the nucleation of AdS^ black holes at Thp = 3/27r. The authors 
of fl^ lay out the ways in which the phase diagram of the weakly coupled theory can 
be matched onto these strong coupling predictions. Various qualitatively different phase 
diagrams can in part be distinguished by the sign of a particular coefficient in the effective 
action for the Polyakov loop U, which is the order parameter for the phase transition. 
The computation of this coefficient requires a three-loop calculation in thermal Yang-Mills 
theory on which is in progress [^ |. 

One of the motivations for the present work was the desire to provide an independent 



check of some pieces of the calculation of [|3^ from a completely orthogonal starting point. 
The one-loop calculation in this paper is equivalent to a two-loop calculation in thermal 



Yang-Mills theory and clearly has some, but not complete, overlap with the work of ||39|| . 
In one sense our calculation contains less information than the effective action for the 
Polyakov loop because we integrate out all degrees of freedom in the theory, including U. 
Moreover, our method is only applicable for temperatures below the Hagedorn transition. 
On the other hand, our calculation contains more information since we have also separately 
calculated, in ( |6.1|) , the one- loop correction to the partition function in the single-trace 
sector. It requires extra work because the 'trace' basis it not particularly natural at 
finite temperature, and the result is slightly messy because of the appearance of 'number 
theoretic' quantities, such as the Euler function (p{n) or the condition k : {k,L) = 1 in 
( ^.35| ). These always disappear at the end of the day in any formula (such as ( ^.141 ) or 
( |6.7|) ) which includes arbitrary-trace operators. 

The finer information present in the single-trace result is deeply related to re- 
cent studies of integrability in A/" = 4 SYM theory because the one-loop correction to 
the partition function is essentially just the trace of the one-loop dilatation operator, and 
therefore encodes the sum of all anomalous dimensions in the theory (sorted according to 
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bare dimension by the powers ofx^o). The one essential subtlety is that only cyclically 
invariant spin configurations correspond to gauge theory operators, because traces of ele- 
mentary Yang-Mills fields are automatically cyclically invariant. Although the spin chain 
Hamiltonian only acts on two neighboring sites at a time, the projection onto cyclically 
invariant states induces an effective long-range interaction on the spin chain which is irrel- 
evant at temperatures near the Hagedorn transition but dominates the partition function, 
and significantly complicates the calculation thereof, at low temperatures. Interestingly, 
the fact that the P5'L(4|4) Hamiltonian is integrable played absolutely no apparent role 
in our calculation. It would be very interesting to understand our result in the context of 
integrability. 

It would also be interesting to consider higher- loop corrections to our results. Although 
the dilatation operator of the full A/" = 4 theory is only known to one loop, in the SU{2) 
sector its precise form is known up to three loops (at the planar level). Furthermore, 
depending on the assumptions (such as integrability) that one is willing to make, one can 
go all the way to five loops 0. 
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Appendix A. The Module Vi 

We tabulate here the SL{4) x SL{2) x SL{2) decomposition of the primary states in 

i - 

the module Vi, using the notation of |^] (where it is referred to as B^-^^ 

2 [l,0,l]{o,o) 



5 [1.0-01(1,0) [0-0,l](o,i) 

2 [0,1,11(1,0) [l'l'0](„_i, 

/ \ / \ 

o [0,l,0](o,o) [0,0,0], [0,2,0](iij [0,l,0](o,o) 

[0,0,2], [0,l,0](i,o) [1,0,1], 1 1, [2,0,0], [O,l,0](o,i) 

/ \ / \ / \ 

/\/ \ / \ / \ 

A mnni [mm [0,0,0], [l,0,l](i,i),[0,2,0](o,o) , „ 

4 [0,0,0](i,o) [0.1.0](i,3) -[0,0,0], [l,0,l](oo) [°'^'°](i,i) [0,0,0](o,i) 

\ / \ / 

9 [0.0,l](i,3) [l'0,0](3,i) 

2 -[0,0,1], [l,l,0](Q,ij -[1,0,0], [O,l,l](i_o^ 



[0,0,0], 3 3^ 

-[o,o,o],[i;o"it 

•■2 ' 2 



,1) 

\ 



^ -[0,0,l](i,i) -[1,0,0](, 1) 



i -[0,0,0](i,i) 

With the help of the SL{A) x SL{2) x SL{2) dimension formula 

dim[/c,p,Q](,^,,,) = -L(A;+p + Q + 3)(A;+p + 2)(p + Q + 2)(A; + l)(p + l)((/ + l)(2ii + l)(2i2 + l) 

(A.l) 

we can immediately read off the partition function for primary states in Vi, 

Vi{x) = Ihx^ + 96x^/2 + 252x^ + 336x^/2 + 210x^ + Ox^/^ - 84x^ - 48x^^/2 - 9x^ (A.2) 
in agreement with the result written in (|5.7| ). 
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Appendix B. Some Details 

In this appendix we sliow liow to obtain the formula ( |5.24| ) from the matrix elements 
of D2 given in 



B.l. The GL{A\4) oscillator basis 

The GL(4\4) oscillator basis for A is realized by a set of four bosonic oscillators a", 
(a, a G {1, 2}) and four fermionic oscillators c" (a G {1, 2, 3, 4}) with the usual relations 

[a°, a^,] = 5^, [b", bp = 5^, {c«, 4} = 51 (B.l) 

and a vacuum |0) annihilated by all of the lowering operators. The only constraint on 
physical states is that they should be annihilated by the central charge 

C=l- ^at^a- + Ib^b- - Ictc^ (B.2) 

The tree-level dilatation operator corresponds to 

Do = l + ^a^a" + hih\ (B.3) 

To consider two letters ^ x ^ we simply have two copies of the above algebra, indexed 
by a subscript (i) G {(1), (2)}. A general state in A(^i) will be labeled by its oscillator 
occupation numbers (a^^j^, a^^^, b^-y 6^^^, c^-s^, c^-y c^-y c^^^) = A"(j-). In this basis we have 

C(^,) = 1 - i(a;,) + al^) + + by - + 4) + 4) + 4)), 

I)o(^) = 1 + i(a;,) + aj,)) + + 4)). 

Then to calculate a trace T^A^i) literally means that we perform a sum of the form 

00 1 

-^(^^ «Jo '''(0 '''W =° '^(0 '4) '^(0 =0 

over all possible oscillator numbers, subject to the physical state constraint. As a check, 
it is straightforward to confirm using this formula and (|B.4|) that 



in agreement with (|2.4| ). 
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Now we consider the action of the one-loop dilatation operator -D2(i2) two sites, 
following ||3^. If we let A| = (aj^, bj^, cj^) schematically denote all of the raising operators, 



then a general state in ^ x ^ can be written as 

.;{/,}) = At^(^^).--At^(^^)|0), (B.7) 

where Si G {1,2} indicates on which site the oscillator acts. The dilatation operator does 
not change the type or number of elementary GL(4|4) oscillators, but can only cause them 
to hop from site 1 to 2 or vice versa according to the rule given in [^: 

^2(12)|S1, • ■ • , Sn; {/i}) = ^ C(n, ni2,7l2l)5(C(i))5(C(2))|si,...,S^;{/J), (B.8) 

<,...,< = 1,2 

where ni2, n2i count the number of oscillators hopping from site 1 to 2 or vice versa and the 
coefficients c(n, ni2, ^21) are given in ( p.26| ). In what follows we will consider a generalized 



operator of the form (|B.8| ) , 

Q|si, . . . , s^; {U}) = Wi, ■ • • , 4; {h}), (B.9) 

with matrix elements q^qi2^q2i^ instead of c(n, ni2, n2i). 

B.2. The combinatorics of hopping 

We are therefore interested in the studying the combinatorics of oscillators hopping 
between two sites. Consider first a toy system with just a single bosonic oscillator a. The 
most general state in ^(1) x ^(2) would then be 

n 

|a(i),a(2)) = (af,))«w(af2))"(^)|0) = n«l.)|0)' ^^'^^^ 

with 

n = a(i)+a(2), = {1, . . . , 1, 2, . . . , 2}. (B.ll) 

a(i) a{2) 

A quantity of interest is the matrix element 



^a(i),a(2)(?'?12,92l) = (a(l) , ^(2) | 1 0(1) , 0(2) ) 

2 n 

= E ?"?r2"92T(a(i),a(2)|Pn«K)|0)' 



(B.12) 
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where Q is defined in ( [B.9| ). The function h counts the number of ways that the initial state 
'^(2)) is mapped to itself, up to a permutation P, under the action of an operator of 
the form ( p.8| ), weighted according to the powers of qi2 and q2i which tell us the number 
of oscillators that have hopped from 1 to 2 or vice versa. An elementary combinatoric 
analysis reveals that 



(5, 512,^21) = ^ {^^^) ("f ) ^'^^^^ ^^^^^ "'^52?' ' 



(B.13) 



= (g?12)"(^U992l)"'^'F(-a(i), -a(2), 1, (?12g2l)"') 



where F is the hypergeometric function. 

A similar analysis can be done for the case of a single fermionic oscillator c. For a 
general state |c(i), C(2)) we find 

^c(,,,c(,,(g,g(i2),?(2i) = {qqi2r''Kqq2ir'^ (l - ^^^"l . (B.14) 

V ?12?21 / 

Since fermionic oscillators only have occupation numbers which are or 1, this formula 
encodes the four cases 

fi'0,0 = 1, ^0,1 = 9921, ^1,0 = g<7i2, fl'1,1 = g^(gi2g2i - i)- (B.15) 

For (7o,o there are no oscillators, so there is no hopping possible. For (70,1 we start with one 
oscillator on site 2, which moves to site 1 giving a factor of q2i- In the final case, gi^i we 
have one fermionic oscillator on each site. They can either stay where they are, or they 
can flip, accounting for the two terms in gi^i. 

For a system with multiple oscillators we simply multiply together the appropriate 
partition functions ( [B.13| ) and ( ]B.14] ) for the individual oscillators. For G'L(4|4) this gives 

^iV(i),Ar(2)(9,9l2,?2l) = {N{l),N^2)\PQ\N[i),N^2)) 

AAA (B.16) 

= n «f2)) n n ^^kd^ ^(2))- 

a=l d=l a=l 

The quantity 

R{w, y- q, gi2, q2i) = TrAxA[Pw''^'''y''''''Q] 

= E 5(C^(l))5(C^(2))«^''°<^'2/''"<^'^iV(,),iV(,)(?,gi2,92l) ^^-1^) 
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then represents a trace over A x A which counts all of the possible hoppings between 
a state |A'"(i), A'"(2)) and its permutation |A'"(2), A'"(i)), weighted appropriately by w to the 
power of the dimension of site 1, times y to the power of the dimension of site 2, times q 
to the power of the total number of oscillators on both sites, times qi2 to the power of the 
number of oscillators hopping from site 1 to 2, times q2i to the power of the number of 
oscillators hopping from site 2 to 1. Concretely, we obtain from ( |B.5| ), ( [B.13| ) and ( [B.14| ) 
the formula 

oo 1 

R{w,y;q,qi2,q2i) = ^ 

"(l)'"(l)'''(l)'''(l)'"(2)'"(2)'''(2)'''(2)-^ (1)' (1)' (1)' (1)' {2)''-{2)''-{2)'^(2) 



S{C^,))6{C^2)){qqi2r^'^{qq2ir^'^w''oMy''o(-^ J] ^ 



a=l 

2 2 

X 

a=l d=l 



4 / a a 

^(1)^(2) 

?12?21 



n ^(-a^i)' -«r2)' 1' (Ml)-') n ^(-^fi)' -^^2)' 1' (Ml)-') 

(B.18) 



where 



a=l a=l a=l 

denotes the total number of oscillators on site i. 

The first step in simplifying (|B.18|) is to rearrange the bosonic sums according to 

E /(4)'4))= E f^h^)^H^)-h^))■ (B-20) 

(For the b oscillators we will use s and t as the new summation variables.) After this 
substitution, the t variables only appear in the arguments of the hypergeometric functions. 
The sum over and t(2) can be done with the help of the identity 

E E ^(~^(1)'~^(2)'1'^)^("^(1) +^(1)'~^(2) +^(2).l'^) 
t{i)=0t(2)=0 (B-21) 

= (1 + S(i))(l + S(2))F(-S(i), -S(2), 2, z). 

The sum over fermionic occupation numbers can be similarly simplified with the for- 
mula 

E /fE^(i)'E^(2))n(i-^(i)^M 

a=l / a=l 

(B.22) 



■^(1) ''^(l) ''^(l) ''^(1) ''^(2) ''^(2) ''=(2) ''=(2) 



F(i)=0F(2)=0 i=0 
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The results (|E:20|) and ( 1^:2^ ) allow (|KT8|) to be written as 



S(l),S(2),i(l),S(2)=0F(i),F(2),J=0 ^-^^ ^ ^ 

X ^l+S(i)/2+S(i)/2^1+S(2)/2+S(2)/2 

X {qqi2r^'^ (??2i)"(^'i^(-S(i), -S(2), 2, -S(2), 2, 

2 

X 

i=l 



(B.23) 

where 

?12?21 

Of course, we are not particularly interested in the quantity R{w,y]q,qi2,q2i)- It 
is simply an auxiliary quantity which counts the possible ways for the operators to hop 
between sites. We are interested in the trace of ( p.8|) instead of the trace of ( p.9|) , which 
we can calculate via the replacement 

{PD2{w,y)) = Riw,y;q,qi2,q2i)\q^q^i2g^2i^^^^^^^^^^^^^ , (B.25) 

where the notation means simply that we expand ( p.23| ) in powers of q, qi2 and §21 and 
then make the substitution indicated for the various powers. 

In order to proceed, we expose the powers of qi2 and q2i in ( [B.23[) by using the identity 



^(-■S(l), -■5(2), 2, -S(2), 2, z) 

^ t S(i)ls(2)l (B.26) 



where JF is the regularized hypergeometric function. Combining ([B.25|) , (|B.26|) and ( [B.23| ) 



finally leads to the formula ( p.24| ) for the desired trace, after some simplification of the 
notation. 
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